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Abstract — The flux upper limits of the diffuse gamma rays, 
from the inner and outer Galactic planes, are revised by factors 
of 4.0^3.7 for mode energies 3^10 TeV, respectively, by using 
the simulation results of the effective area ratios for gamma-ray 
induced showers and cosmic-ray induced ones in the Tibet air 
shower array. In our previous work, (Amenomori et al., ApJ, 
580, 887, 2002) the flux upper limits were deduced only from 
the flux ratio of air showers generated by gamma rays versus 
cosmic rays. The details of the simulation are given in the paper 
(Amenomori et al.. Advances in Space Research, 37, 1932, 2006). 
The present result using the same data as in ApJ suggests that 
the spectral index of source electrons is steeper than 2.2 and 2.1 
for the inner and outer Galactic planes, respectively. 

I. Introduction 

Diffuse gamma rays in MeV^GeV energy region from the 
inner Galactic (IG) and outer Galactic (OG) planes observed 
by EGRET (Hunter et al. 1997) [1] show a sharp ridge both 
along the IG and OG planes. The EGRET flux is about 3 times 
higher than COS B data (Mayer-Hasselwander et al. 1982) 
[2] in several GeV, although the flux is consistent with the 
conventional calculation (Dermer 86) [3] in £^ < 1 GeV. The 
EGRET excess above 1 GeV has been tried to explain by 
some models; a hard source electron spectrum of index (3=2.0 
by Pohl & Esposito (1998) [4], hard proton spectra by Mori 
(1996) [5] and Webber (1999) [6], and an additional secondary 



electrons and positrons raising from the cosmic-ray collisions 
with ISM by Strong et al. (2004) [7]. 

In higher energy region theoretical calculations have been 
given by Porter & Protheroe (1997) [8] and Tateyama & 
Nishimura (2001) [9] for the inverse Compton (IC) gamma 
rays, and by Berezinsky et al. (1993) [10] for tt^ 27 
process through cosmic-ray interaction with ISM. The most 
experimental data in higher energy region gave only flux upper 
limits except the definite flux by Milagro [11] at 3.5 TeV 
for IG. Not only the absolute flux but also the flux upper 
limit are both important for restriction of theoretical models. 
In this paper, using of the detection area ratio of the Tibet 
array between gamma rays and galactic cosmic rays in the 
simulation, we revised and decrease the flux upper limits in 
our previous paper (1002) [12], and compare with the Milagro 
result. The details of the simulation is described in our recent 
paper (Amenomori et al. 2006) [13]. 

II. Simulation of effective areas 

Shower size of primary gamma-ray induced showers is 
about three times larger than galactic cosmic-ray induced ones 
in average at the d epth of 606 gm~^ of the Tibet array for the 
multi-TeV energy region. Hence, the effective area of the array 
is larger for gamma rays than cosmic rays. Figure 1 shows the 
layout of the Tibet III air shower array. The Tibet II array is 




Fig. 3. Assumed primary spectra and distributions of triggered events in the 
simulation for IG plane. 
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Fig. 1. Layout of the Tibet III array in Yangbajing at the stage in 1999 
2001. 
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Fig. 4. Simulation results of effective areas of the Tibet III from IG plane 
for primary gamma rays, protons and all cosmic rays. 



Fig. 2. Tibet III exposure map in the galactic coordinates for the zenith 
angle with ^ < 50°. 



the one excluding detectors with the open squares. Figure 2 
is the exposure map of the Tibet III array for the zenith angle 
of < 50°. 

Figure 3 shows differential energy spectra of triggered 
gamma induced showers and cosmic-ray induced ones, assum- 
ing both spectral index of 2.6 in the simulation. We can see 
the mode energy of triggered gamma rays is 1.5^2.0 times 
smaller than cosmic rays for the same trigger condition, and 
the effective area ratio to be about 7 in multi-TeV region. 

Figure 4 shows the average advantage factor of the effective 
area for gamma rays versus cosmic rays is 4.0 for E^mode — 
3 TeV and 3.6 for 10 TeV in average of gamma-ray spectral 
indices of P = 2.2 ~ 2.8. In this paper the mediate value 
P=2.5 is employed because of its weak dependence on the 
spectral index. 

The significance of the excess sigma (a), in the Table 1, of 
TeV gamma rays from the IG and OG planes implies a simple 



formula of {E — B)/ V^, where E is the number of events 
on-plane and B is the background number of events estimated 
from neighboring bins around the on-plane. Those evaluated 
values in the previous paper [12] are tabulated together with 
the simulated effective area ratio of gamma rays vs. cosmic 
rays and the revised upper limits in a case of small change of 
source electron spectral index of (3 from 2.4 to 2.5. 

III. Results and Discussions 

The present simulation, giving a larger effective area of the 
Tibet array for gamma rays than the galactic cosmic rays, 
decreases the flux upper limits of diffuse gamma rays by a 
significant factor, as given in Table 1 . The original data in the 
previous paper (Amenomori et al. 2002) [12] were obtained 
by the Tibet III array with inner area of 22,050m^ at 3 TeV, 
and by the Tibet II array with 28,350m2 at 10 TeV. The details 
of the simulation of the effective area ratio between gamma- 
ray /cosmic-ray initiated air showers are described for the Tibet 
III array, in the recent paper (Amenomori et al. 2006) [13]. In 
Fig. 3 is shown the detected shower event distribution in the 
simulation, and in Fig. 4 is also shown the effective area for 



TABLE 1 Effective area ratio for gamma rays vs. cosmic rays and gamma-ray flux upper limits. 



Galactic latitude range 


|6| < 2° 


|6| < 5° 


Air shower array 


Tibet III 


Tibet II 


Tibet III 


Tibet II 


-^mode 


3 TeV 


10 TeV 


3 TeV 


10 TeV 


Inner or Outer Galactic planes* 


IG 


OG 


IG 


OG 


IG 


OG 


IG 


OG 


vjlgjllllll^cllll^C Ui CAl^CSS \(J ) 


+2.52 


+0.25 


+1.71 


-0.63 


+1.88 


+1.78 


+0.81 


-0.66 


Flux ratio of 7 rays vs. cosmic rays"^ 


















/^(1(7)//CR = I/a/B (xlO-4) 


1.95 


1.16 


2.43 


1.45 


1.23 


0.737 


1.54 


0.936 


99%CL revised flux upper limits 


















J(> E^) (10-iicm-2s-isr-i) 


57.3 


19.4 


7.81 


2.50 


31.4 


18.4 


3.95 


1.60 


99%CL upper limits {(3=2A)^ 
E'^dJ{> E^)/dE^ 


9.6 


3.3 


4.0 


1.3 


5.3 


3.1 


2.0 


0.83 


(xlQ-^cm-^s-isr-^MeV) 


















Effective area ratio of 7-ray/CR 

'Seff(7)/'^eff(CR) 


4.0 


3.7 


4.0 


3.7 


99%CL revised upper limits (^5=2.5) 


















E^dJ{> E^)/dE^ 
(xlQ-^cm-^s-isr-^MeV) 


2.6 
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t are referred to Amenomori et al. ApJ., 580, 887 (2002) 

* Ranges of IG are 20° < / < 55°, and of OG 140° < I < 225°. 



primary gamma rays, protons and all cosmic rays. 

Figures 5 and 6 show the revised flux upper limits, for IG 
and OG planes, at £^mode=3 TeV (T3: for the Tibet III array) 
and at 10 TeV (T2: for the Tibet II array). In these figures the 
EGRET data (Hunter et al. 1997) [1] of the Galactic latitude 
width of \b\ < 2° are shown, and also shown the upper limits 
by Whipple (W) (LeBohec et al. 2000) [14] with 99.9% C.L. 
and HEGRA (H) (Aharonian et al. 2001) [15] with 99% C.L., 
though both at a small sky region around the galactic longitude 
of / = 40°, and HEGRA-AIROBICC (Ha) (Aharonian et al. 
2002) [16] and CASA-MIA (C-A) (Borione et al. 1998) [17] 
both with 90% C.L.. Theoretical curves of inverse Compton 
are given by Porter and Protheroe (1997)[8] (PP2.0 and PP2.4 
in figures), and by Tateyama and Nishimura (2003) [9] (TN 2.0 
and TN2.4), where 2.0 and 2.4 are assumed source electron 
spectral indices. Theoretical curves arising from tt^ 27 
decay are also given by Berezinsky et al. (1993)[10] (BGHS). 

When the observed gamma-ray spectra with the spectral 
index of 2.5 is adopted in this paper, the revised results can 
give a strong suggestion that the spectral indices of source 
electrons for the inverse Compton (IC) are steeper than 2.2 in 
the IG plane and also 2.1 in the OG plane in comparison with 
the theoretical calculations of IC. 

Recently, Milagro (M) obtained the definite flux at 3.5 TeV 
from IG plane and the upper limit from OG plane. Both are 
a few times lower than ours. In their IG data are involved the 
range of 40° < / < 100° and \b\ < 5°, though our longitudinal 
IG range of 20° < / < 55°. Their result suggests that the 
diffuse gamma rays from the Galactic plane can be interpreted 
not only by IC but also by tt^ 27 below 10 TeV. The 
longitudinal range of Milagro is larger than our range of 20° < 
/ < 55° and the number of events is 3.6 times larger than ours 
for IG with \b\ < 5°. The advantage factor of hadron rejection 
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Fig. 5. Diffuse gamma rays from the IG plane of 20° < / < 55° for the 
belt with |6| < 2°. 
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Fig. 6. Diffuse gamma rays from the OG plane of 140° < I < 225° for 
the belt with |6| < 2°. 
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Fig. 7 Integral flux upper limits of Tibet data with |6| < 5° plotted 
on the figure referred from the Milagro's paper PRL, 95, 251103(2005). 



of Milagro is estimated to be 0.45/a/oT ~ 1.45. Thus 
the effective area ratio of detection of gamma-ray showers 
vs. hadron showers is 1.45>/3^ ~ 3 times better than our 
99%C.L. upper Hmit corresponding to 4 ~ 5cr. Their prots of 
Milagro data are lower by a factor of about 4 than ours. This 
might be due to an odd analysis method of "time swapping". 
Otherwise the effective area estimated by Milagro is large. 

It seems that the recent Milagro result [11] rather claims 
TT^ 27 process by cosmic-ray particles according to the 
interpretation that their data is on the straight line extrapolated 
from the data above 10 GeV of EGRET (Hunter et al. 
1997)[1]. Figure 7 is referred from Milagro paper (Atkins et 
al. 2005) [11] and includes our present revised upper limits for 
\b\ < 5°. This range of the Galactic plane covers the Cygnus 
region in which a lot of complicated point sources or diffuse 
gamma-ray sources seem to be involved. This is the reason 
why we have excluded this region from the analysis of the 
inner Galactic plane, and we have been examining the Cygnus 
region separately from IG plane. 

Recently, we reported our anisotropy result for the multi- 
TeV galactic cosmic rays in the northern sky [18]. In that 
paper, a broad hill excess of cosmic-ray intensity is found in 
the Cygnus region in 13.3(7, and clear excesses are observable 
for different energy ranges in the same region. 

IV. Conclusions 

The Tibet AS7 data collected between 1997 and 2001 
have been reanalyzed with improved background efficiency 
calculation by about a factor of 4. Assuming an index of 
2.5 for gamma-ray spectrum when calculating gamma-ray 
efficiency, 99%C.L. flux limits are derived which suggests 
that the corresponding source electron spectrum for the inverse 
Compton (IC) is steeper than 2.2 for the IG plane and 2.1 



for the OG plane. However, those gamma-ray flux limits 
are too high to constrain proton source models. In the near 
future, the broad hill of excess in the Cygnus region will be 
further studied with more available statistic. As the current 
Tibet air shower array is not able to discriminate the gamma- 
ray induced air shower from the hadron-induced one, we are 
planning to add dozens of large water cherenkov detectors in 
the Tibet array, to reject the hadron initiated air showers by 
using the information of secondary muon. Our gamma ray 
sensitivity is foreseen to be improved by a factor of 4^10 
at 10^100 TeV, while at PeV energy we will be free from 
background [19]. 
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